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Abstract 

Integral bases, a minimal set of solutions to Ax < b, x € Z™ that generate any other solution 
to Ax < b, x £ Z n , as a nonnegative integer linear combination, are always finite and are at 
the core of the Integral Basis Method introduced by Haus, Koppe and Weismantel. 

In this paper we present one generalization of the notion of integral bases to the nonlinear 
situation with the intention of creating an integral basis method also for nonlinear integer 
programming. 

1 Introduction 

In the past fifty years many efforts have been undertaken to study linear integer optimization prob- 
lems from different mathematical and algorithmic viewpoints. As a result, a basic understanding 
of the geometry of integer programming problems defined by linear equations and/or linear in- 
equalities is present today. This knowledge has been partly turned into algorithmic tools to tackle 
discrete optimization problems in practice. 

The attempts to study the geometry of integer points in polyhedral sets are based on two basic 
mathematical concepts. One is the notion of a lattice. More precisely, a basis of a lattice L is a subset 
of linearly independent vectors that allows one to generate all points in the lattice with respect to 
taking integer linear combinations. The geometric properties of particular bases in lattices made 
it possible to design algorithms for solving specific linear integer programming problems, mainly 
problems without lower and upper bounds on the variables and linear problems with a fixed number 
of discrete variables Ej • The notion of bases of a lattice can be further refined so as to yield 
so-called integral generating sets for cones and polyhedra. Roughly speaking, integral generating 
sets extend - besides lattices - the notions of extreme points and rays in polyhedra and cones 
to integer points in such sets. More precisely, an integral generating set for a set S C Z" is a 
subset of S with the property that every member of S can be represented as a nonnegative integer 
combination of the elements in S. Of course, S itself constitutes an integral generating set of itself. 
The key question is to detect an integral generating set that is finite and minimal with respect to 
inclusion. This immediately raises the question to characterize those sets S of lattice points that 
possess a finite integral generating set. This question is answered in Section |3] of this paper. 

Indeed, integral generating sets have important implications for the theory of linear integer pro- 
gramming. Most importantly, optimality conditions for integer optimization problems can be de- 
rived through integral generating sets. Such sets also provide a basic understanding of integral 
polyhedra and totally dual integral systems of inequalities 2 . Last but not least they play a cen- 
tral role in the development of integer simplex type methods of linear integer programs |4|. In fact, 
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it is quite obvious to see that if a finite integral generating set for a discrete set of points is avail- 
able, then we can reformulate the problem of detecting a particular element in S as the problem 
of detecting a nonnegative integer multiplier associated with the new representation through an 
integral generating set. Integral generating sets therefore allow a new representation of the same 
set S in some other space. The beautiful fact is that if we start off with a set S that is the feasi- 
ble region of an integer linear program in nonnegative variables, then also after reformulation the 
new optimization problem happens to be a linear integer program in nonnegative variables. This 
follows simply from the fact that the integral generating set enables us to express every point as a 
nonnegative integer combination. 

Suppose now that S does not have a finite integral generating set. This in fact may happen even 
though S corresponds to all the integer points in a polyhedron. In particular, if the constraints 
defining S are not linear, even in quite restrictive cases S does not possess a finite integral generating 
set. 

Then the idea to use integral generating sets for reformulation issues is not possible, because the 
generating set is infinite. We can simply not write down any finite representation of the reformulated 
problem. In order to cope with this scenario, it requires to generalize the notion of an integral 
generating set from the linear case to a nonlinear setting. We will refer to such sets as integral 
function bases, since they enable us to derive representations by means of nonnegative polynomial 
combinations instead of nonnegative linear combinations. This is the central topic of Section 
In turn, our generalization allows us to formulate optimality conditions for integer polynomial 
programming problems. 

We also analyze the situation when an integral function basis for the integer points in a complicated 
semi-algebraic set is replaced by the condition of being members of a relaxation of the semi- 
algebraic set itself. This question is in particular motivated by the design of pivoting type methods 
for polynomial integer and mixed integer programming, a topic that we regard as theoretically and 
practically challenging, but important. 

2 Integral Bases 

Let us start by defining the notion of an integral basis. 

Definition 2.1 Let S C U 1 . Then we call T C S an integral basis of S, if for every s E S there 
exists a finite (integer) linear combination s = Y)ajti with ti £ T and oti G Z+. 





Figure 1: Minimal integral bases of two sets of lattice points 
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Note that an integral basis of S is allowed to contain elements only from S itself! With this 
definition, Bcrtsimas and Weismantel £Q showed the following characterization of which rational 
polyhedra (or more precisely the integer points in such polyhedra) have a finite integral basis. 



Theorem 2.2 (Bertsimas and Weismantel |1J) For A £ Z dx " and b £ 7L d , define the sets 
P = {x € M.% : Ax < b}, S = P n Z", and C = {x e : Ax < 0}. 



(a) There exists a finite integral generating set of S if and only if S contains all but finitely many 
integer points in C H U\ . 

(b) If a finite integral generating set of S exists, then there is a unique integral basis of S . 

Note that for b = 0, this theorem simply states existence (and uniqueness) of a Hilbert basis for 
the (pointed) rational polyhedral cone {x e K7? : Ax < 0}. 

Now let us give a novel and more general characterization of which sets of lattice points have a 
finite integral basis. As we do not make any structural assumption on the set of lattice points, 
we have to be cautious to check whether the integral bases that we construct do indeed consist of 
lattice points from our original sets only. 

Theorem 2.3 Let S C Z™ be any set of lattice points in Z™. 

(a) S has a finite integral basis if and only if C — cone(S') is a rational polyhedral cone. 

(b) If the cone C = cone(S') is rational and pointed, there is a unique finite integral basis that is 
minimal with respect to set inclusion. 

Proof. Let us start showing part (a). If G = cone(S') is not a rational polyhedral cone, S cannot 
have a finite integral basis G C S, since C = cone(S') = cone(G) would be a rational cone, 
contradicting our initial assumption on C. 

Now we show the remaining claim that S has a finite integral basis if C — cone(5) is rational by 
explicitly constructing such a finite basis. It should be noted that this integral basis need not be 
minimal. 

First, let us triangulate C into (finitely many!) simplicial cones C\, . . . , Note that we can and 
do choose such a triangulation for which the generators of the cones Cj are also among the (finitely 
many) generators of C. Thus, as C, each cone Ci is generated by (finitely many) elements Si of S. 
It remains to show that for each rational simplicial cone Ci = cone(S'i), the set d PI S has a finite 
integral basis G;. Then the union of all G,, i = 1, . . . , k, is clearly a finite integral basis for S. 

For d — cone(Si) and Si = {v\, . . . , v r }, consider the parallelepiped 



As Fi is bounded, F contains only finitely many lattice points {/i, . . . , ft} in Z n . Moreover, Gj PlZ 
is the disjoint union of the following t sets F% , . . . , F t with 
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We construct now a finite integral basis for Cj n S. 

Consider any Fj, j — 1, . . . ,t. As Cj is a simplicial cone, each point in Fj has a unique representation 



To construct a finite integral basis for Fj n S*, consider the set 4>j{Fj Pi 5) C Zj_. By the Gordan- 
Dickson Lemma, there are only finitely many points . . . , g p } that are minimal with respect to 
the partial ordering < defined on 11^. Thus, each point 4>j(Fj n S) can be written as a positive 
integer linear combination of gi, . . . , g p and of the unit vectors e%, . . . , e r . Thus, every element in 
Fj (~l S is a positive integer linear combination of _1 (gi), . . . ,<j>~ 1 (g p ) 6 Fj H S together with 
_1 (ei), . . . ,</) _1 (e r ) £ C 5. Let Gjj denote the set of all these vectors. Clearly, the union Gi 
over all Gij, j = 1, . . . ,t, forms a finite integral basis for Ci f) S, and claim (a) is proved. 

Let us prove claim (b) now. As cone(S') is pointed, there is some vector c 6 K n such that {x 6 
R™ : c J x < 0} n cone(S') = {0}. Assume that U = {ui, . . . , u r } and V = {vi, . . . , v t } are two 
different inclusion minimal integral bases of S. Moreover, assume that w.l.o.g. u\ ^ V . Minimality 
of U implies that u\ cannot be written as a positive integer linear combination of elements in 
U \ {ui}. However, as V is an integral basis of 5* and u\ £ S, there is a nonnegative integer linear 
combination u\ = Y)j—\ ctjVj. Clearly, as u\ ^ V and as the coefficients are nonnegative integers, 
we have c J Vj < c J ui whenever a.j > 0. As also U is an integral basis of S and as all Vj € S, there 
are nonnegative integer linear combinations Vj = flijiii. Moreover, c J Ui < c T Vj whenever 

0i t j > 0. Plugging these representations into u\ = 0L j v ji we a representation of u\ as a 

nonnegative integer linear combination of elements in U . However, by construction, they all have 
a scalar product with c that is strictly less than c J U\. Thus, we have written u\ as a nonnegative 
integer linear combination of elements in U \ {ui}, a contradiction to our assumption that U is a 
set inclusion minimal integral basis, and the claim is proved. □ 

Note that for sets of the form (x £ Z" : Ax < 0}, Theorem 12.31 again simply states existence of 
finite Hilbert bases for rational polyhedral cones and uniqueness of the minimal Hilbert basis if 
the cone is pointed. It is easy to show that the minimal Hilbert basis of a cone must consist of 
lattice points from the fundamental parallelepiped and is thus finite. The tricky part for the proof 
of Theorem 12.31 was the fact, that not all points of this parallelepiped could be assumed to belong 
to S. The two examples in Figure ^ page [3 already illustrate this difficulty. 

Let us now show how Theorem 12.31 implies the special case, Theorem 12. 21 

Proof of Theorem 12.21 Let us show part (a) first. If S is finite, nothing is left to show. Thus, 
assume that S is not finite and therefore also C ^ {0}. Assume that S contains all but finitely 
many integer points in C (~1 Z!f . In particular, S contains an (integer) point of every extreme 
ray of C. By Minkowski's theorem, we have conv(P nZ") = conv(G) + C, where G C S is the 
set of extreme points in conv(P n 1 n ) C P. (Since P does not contain a line, G ^ 0.) Thus, 
cone(S') = conv(G) + C, as G C S and as S contains an (integer) point of every extreme ray of C. 
Consequently, cone(5) is a rational cone and thus has a finite integral basis by Theorem l2.3l 

Now assume that there are infinitely many integer points in C that do not belong to S. In particular, 
^ P as otherwise = AO < b implying C C P and thus C f) Z" C S. Assume for the moment 
that each extreme ray of C contains a (nonzero!) point of P. Fix any extreme ray of C and let 



as fj + E r r . 
by 



ctjVj implying that there is a one-to-one correspondences 4>j between Fj and TL given 
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v £ C n P be a point on this ray. Then any point w — kv, k > 1, on this ray must belong to P. 
This follows from Aw = k(Av) = (k - l)v + v < (k — 1) • + b = b and Av < 0, Av < b, and k > 1. 
Therefore, we Pas claimed. By convexity of P, P must contain the convex hull H of all these 
half- lines {kv : k > l,k GM}. As C\H is bounded, only a finite number of integer points in C can 
lie in C \ H . As S = P D Z n , this implies that only finitely many integer points C can lie outside 
of S, contradicting our initial assumption on C. This implies that there must be an extreme ray 
R of C that does not contain any point of P. 

We now show that cone(S') cannot be a rational cone, and the result follows again by Theorem l2.3l 
Assume on the contrary that cone(S') is a rational cone. By convexity of S 1 , every ray in cone(S') 
has a nontrivial intersection with S and thus also with P. This implies that the extreme ray R 
does not belong to cone(5). As cone(S') is rational, there exists a finite (rational) description 

cone(S') = {x £ R n : cjx < 0,i = 1, . . . ,p}. 

Let v be any rational vector with R = cone(w). Then v cone(S') implies that there is some index 
j such that cjv > 0. Now consider any integer point w £ S. As S — conv(G) + C, all integer 
points on the half-line {w + av : a > 0} belong to S. Moreover, as v is a rational vector, there are 
infinitely many integer points on this half-line. However, as cjv > 0, we have cj(w + av) > for 
sufficiently large a, implying that there are integer points of S that lie outside of cone(S'). This 
contradiction shows that cone(S') is not a rational cone and part (a) is proved. 

As part (b) of our claim follows now immediately from part (b) of Theorem 12 .31 nothing is left to 
show. □ 

A natural question that we may ask is, whether there are other special cases of interesting sets of 
lattice points that have a finite integral basis by Theorem 12. 31 

One natural guess would be the integral points in a convex region. However, convexity alone is not 
enough to ensure that cone(S') is rational, as can be seen by looking at a polyhedral cone with 
irrational generators. Thus, some notion of "rational generators" of the region should be defined. 

With this in mind, we may try to look at sets that are parametrized by convex polynomials that 
have rational coefficients only. Again, there is a simple counter-example. For the lattice points S 
in the parametrized set 

{(;) e - ! (;)-(^) + (?)-' £ 4 

we easily see that cone(S') is not rational, see Figure [3 



Figure 2: cone(S') is not rational. 
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We conclude that even convexity and rationality of generators is generally not enough to ensure 
finiteness of an integral basis. It can be shown that under the assumption that if in addition the 
given set itself is convex and that if it contains the unit vectors of the positive orthant, a finite 
integral basis does exist. 



3 Nonlinear Integral Bases. Definition and Motivation 

Theorem 12.31 characterizes when linear integral bases exist. What can we do if the conditions of 
the theorem do not hold? For instance, if we consider the set S = {(x, y) £ "L\ : y > 1}. 



(0.0) 

Figure 3: Example of an infinite integral basis 

In this case, the set cone(S') is not finitely generated, and thus there does not exist a finite integral 
basis of S. For obtaining a finite representation in this example, it becomes necessary to extend the 
notion of an integral basis to - what we call - an integral function basis. Our goal then becomes 
to identify sets of points that have a finite integral function basis. In the following we consider sets 
S = {y E R™ : y = g(X), A e Z^} flZ™, where g : R d — ► M" is a vector of functions with components 
gi : R d R, i = 1, . . . , n, and with g(Z d ) C Z n . 

Note that when all gi are linear functions, our set S corresponds to the lattice points of a rational 
polyhedral cone. Other possible functions are polynomials in Z[A], certain stair-case functions, or 
even suitable combinations of all 3 types. 

Example 3.1 The function g given by 




defines a semi- algebraic set C — {y € ]R 2 : y = g(X),X £ Z 2 ^}, see Figure^ In Cartesian coordi- 
nates, C can be described by C = {(x, y) £ M 2 : x — y 2 > 0, x, y > 0}. 

Of special interest to us will be the lattice points inside semi-algebraic sets. 

Definition 3.2 Consider a set S C Z". Let sets Ti C Z™ be given where each Ti is described in 
the form Ti := {fi(ti) : ti £ Z!j?} with a polynomial function fi : — > Z". 

Then we call such a family {Ti} an integral function basis of S, if for every s S S there exists 
a finite representation, s = ^2 fi(ti), with ti € Z™ 1 and fi(ti) G S . 



G 



Figure 4: Semi-algebraic set C with its lattice points 



If we allowed only linear functions /, and if S are the lattice points in a rational polyhedral cone, 
this definition coincides with the definition of a Hilbert basis. 

If we reconsider the example with S — {(x,y) G 7!?+ : y > 1}, we see that the following set T\ 
defines an integral function basis of S: 

Ti = {/(A, M ) = (A,l+M):A,MeZ + }. 

Example 13.11 cont. Let us consider again the semi-algebraic set C = {(x, y) £ I 2 : x — y 2 < 

0,x,y> 0}. An integral function basis of S = CC\J? is given by {T 1} T 2 } with 7\ = {(0, t):t€ Z+} 
and T 2 = {((x+l) 2 — s,x+l) : x,s € where the parameters x and s need to satisfy s < (x+1) 2 
to guarantee ((x + l) 2 — s, x + 1) G S, see Figure The only lattice point in C that cannot written 
as a sum of a lattice point in T\ and a lattice point in Ti is the origin. This special point, however, 
can already be represented by T\ alone. 



Figure 5: Integral function basis of a semi-algebraic set 

Note that the constraints that are needed to encode the condition (x+1) 2 — s,x + l) G S have the 
same maximal degree as the original constraints. On the other hand, we can also this condition by 
s < 2x + 1, see FigureEl The latter representation should be preferred, since this description of the 
set S using the integral function basis involves only linear constraints in contrast to the quadratic 
constraint in the description above. □ 

In the following, we outline a fundamental application of integral function bases for nonlinear 
integer optimization problems. It turns out that one can derive an optimality criterion for a linear 
integer program with a polynomial objective function. This criterion is a natural generalization of 
what Graver proved in the fully linear (integer) setting 

Theorem 3.3 Letp(z) be <xny polynomial in . . . , z n 1, and letAeZ dxn and b G Z d . Consider 
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Figure 6: Integral function basis of a semi- algebraic cone 



the (possibly nonlinear) integer optimization problem 

max{p(z) : Az = b,z > 0,z G Z n }. 
Let W l G l Z nxWi , i = 1, . . . 2™, denote the extreme rays of the cones 

cone(VT) = {x G K" n O, : Ax = 0}, 

Thus, every point in this cone can be written as 



where Oi, . . . , G>2" denote the 2" orthants of I 
a linear combination z — W l X, \ > 0. 

Assume that zq is a feasible integer solution to Az 
following vector of nonlinear functions, 



b, z > 0. For each i = 1, . . . , 2™, define the 



'(A) 



/ 51(A) \ 
\ f (A) J 



p(z + W t X)-p(z Q ) 



Let {T[, . . . , Tfi. }, with Tj := {fj(tj) : S R™*' 3 }, &e an integral function basis for the (inte- 



ger points in the) semi-algebraic set C % — {y S 



nn+l 



: V = .° l (A),A > 0}. Define by Sij the 
n e cone(W /l ). (Herein, [fj(tj)]i,..., n shall 



semi-algebraic set that encodes the conditions [fj(tj)]i t . 
denote the vector of the first n components of fj{tj).) 

Then zq is optimal if and only if for every i = 1, . . . , 2™, the following condition holds: 

[/jfe)]n+i < for all t 3 e Z^' 3 n S id with [f){t )] k > -[z ] k , for all k = 1, . . . , n. 

Proof. Assume that there is a better feasible solution z\ that has an objective value p{z\) > p(zo). 
Consider the difference vector v := Z\ — zq, which lies in one of the 2™ orthants Oi of R n . Therefore, 
we are looking for v G Oi with Av = 0, z + v > 0, and p(z n + v) — p(zo) > 0. Clearly, the set 
{z £ Oj : Az — 0} forms a pointed rational cone, generated by the columns of W l . Thus, v — W l X 
for some A G KT* and hence 

W l X 

p(z + W i X)-p(z ) 

is an integer point in the semi-algebraic set C l = {y G M n+1 : y = g z (X), A > 0}. Using the integral 
function basis of this set, there is a representation 

= E/i<y = i;( [ fe 



W l X 

p(z Q + W l X)-p(z ) 
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with t) G Sij n Z n « and /j'(tj) G C\ 

As p(z + VP A) — p(zo) > 0, there must be some j G J* with > 0. We claim that the 

first n components of fj{tj) G C 1 form an improving integer vector for zq , possibly different from 
the vector v that we decomposed. 

As [fj{tj)]n+i > 0; ^e only thing left to show is that the components [fj(tj)]i, ■ . ■ , [fj(tj)] n of /jfe) 
lie above the lower bounds, i.e., [fj{tj)}k > — [zo]k for A; = 1, . . . , n. But this can be seen as follows. 
By construction, cone(VP) G O l , implying [/j(£y)] G O l for all j. Thus, the components of 
z o + [fj(tj)]i,...,n lie between the components of zq and of zq+v — z\, and are therefore nonnegative. 
The converse direction is obviously true. □ 

Clearly, one would wish that searching for an improving vector in each of the Ti is simpler than 
searching for an improving vector in C. 

The set Ti = {A-/i:A,^£ Z™ } always forms an integral function basis for any set S C Z™, where 
2n parameters are needed to describe Ti . The following theorem bounds the number of parameters 
needed in the and thus gives a sufficient condition (together with a construction) of when an 
integral function basis with less parameters in the description of each Ti exists. 

Theorem 3.4 Let SCZ", v\, . . . , Vk G Z n one? Zei C = cone(i>i, . . . , Vk) be a rational polyhedral 
cone with S C C. T/ien 5 /ias an integral function basis in which the appearing sets Ti involve at 
most k + 1 parameters. 

Proof. First observe that 5 C C implies cone(S') C C and therefore cone(S' U . . . , Vk}) = 
C. Thus, by Theorem 12.31 there is a finite integral basis {hi, . . . , h s , vi, . . . , Vk} for the set S U 
{yi, . . . , Vk}- If we set in addition ho = 0, we can see from the proof of Theorem 12.31 that every 
point v G S can be written as v = hi + Y)j—i \jVj for some i G {0, . . . , s} and for some nonnegative 

integers Ay. This last condition in fact states that the sets Tj = {hi + J2j=i ^ij v j '■ e ^+}U{0}, 
i = 0, 1, . . . , s, form an integral function basis for S*. □ 

Remark 3.5 It should be noted that we may strengthen the above theorem if some or all of the 
cone generators Vj lie in S. If Vj G S, then each set Ti = {hi + Y)j—i '■ Xi E Z^} U {0} can 

in fact be decomposed into the sum of T[ — {hi + J2je{i fc}\{jo} ^ i i v i '■ ^ ^ Z^ -1 } U {0} and 

T i = {*ijo v h : X ijo 

Iterating this process for all cone generators vj that lie in S gives a new integral function basis for 
S with fewer parameters appearing in the description of the sets Ti. In fact, if all Vj lie in S , that 
is if C = cone(S') ; the integral function basis for S simplifies to sets Ti that all contain nonnegative 
integer multiples of a single lattice point of S . Thus, we have recovered the statement of Theorem 
\2.<A the existence of a finite integral basis if cone(iS) is rational. 

The following example demonstrates that splitting the set S into finitely many subsets may also 
decrease the maximum number of parameters needed in the description of the Tj's. 

Example 3.6 Consider the set S — {(x,y) G Z 2 : — y < x < y 2 }. As this set is contained in the 
rational cone spanned by e\ and —e\, we conclude by Theorem \3.J\ that S has an integral function 
basis, in which each Ti is described by at most 2 + 1 = 3 parameters. 
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However, if we split the set S as 

S = S' + S" = {(x, y) S Z 2 : -y 2 < x < 0} U {(x, i/)eZ 2 :0<x< y 2 }, 

we see that S' C cone(— e\, e-i) and S" C cone(ei, e^. Since the cone generator ei is an element of 
the sets S' and S" , respectively, we find integral function bases for S' and S" , in which each Ti is 
described by at most 1 + 1 = 2 parameters. Putting both together, we arrive at an integral function 
basis for S with the same property. □ 



As we have seen above, every set of lattice points in Z" admits a representation via an integral 
function basis. Even under the assumption that we have found a nice integral function basis for a 
particular problem instance, that is, one that has only few parameters in the description of the Tj, 
we are faced with a new problem to be solved. 

Suppose we want to maximize a (polynomial) function p(x) over the lattice points in a semi- 
algebraic set C. Knowing an integral function basis {Ti : i 6 /}, we can use the representation 
x = J2iei /•*(*»)> e f° r all x G C n Z™ to rewrite the problem as 

max | p (yi : s /*(**) ecn z "' u e z + } • 

While the condition fi{U) e C n Z™ often follows immediately from fi{t t ) G C nZ n for all 

i E I, these latter conditions involve descriptions by polynomials of the same degree as in the 
description of C and are thus still hard to deal with. Finding ti G Z™* with /(i,) G CD Z™ even 
only for a single i (as needed in Theorem 13. 3|) is as hard as finding a point in C, at least from a 
complexity point of view. 

Thus, an integral function basis with the additional property that Tj C C D Z™ for all i G I would 
be desirable. Then G C n Z" for all i G I and for all ^ G Z™ ; would hold automatically. 

For this, of course, a nonlinear description for the Tj is needed, in contrast to the rather nice and 
simple description guaranteed to exist by Theorem 13.41 

In the following, we relax the condition /j(tj) G Cfl Z" and allow a correction term that may lie 
outside, but which is bounded by polynomials of strictly smaller degree than the given polynomials. 

Theorem 3.7 For every semi- algebraic set 

C := {x G R n :3y>0 with x = g(y)}, g : R d -> K", 
there exists a set of functions 

{g h g u :W l ^ M"} 

with maxdeg((7;), maxdeg((7 u ) < maxdeg(<7) and such that for every point x G C n Z™ i/iere exists a 
A G Z^ and a point v x G Z™ wit/i a; = g(A) + w x and with gi(X) < v x < g u (A). 

Proof. Choose any x G 5 := C n Z™. Then a; = g(t/) for some y G R™. Now define A := [y\ 
component- wise and let v x = x—g(X) and h = y — X. We will now construct functions gi : R n — > M n 
and g u : R™ — > M™, with the desired properties. 
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Let D = maxdcg(g). By multivariate Taylor expansion, we get for j = 1, . . . , n: 

^)= 9 0)(A + / l ) = 5 0)(A) + f:I. £ ^P-^- 

1=1 aeZ" + : 



Therefore, 



ah = i 



i\ ^ dx c 

e e Z™ : 

alh = i 



Note that maxdeg ^7^3-^ < maxdcg(g) and that < h < 1 by construction. 

This sum is a polynomial in A and h, that is, it is a sum of terms c a jiX a h^ . Since all A > we can 
use < hi < 1, for all i, to bound the expression c a ^X a hP by 

< CcpX"^ < c a , p \ a 

if c a ,/3 > and by 

c a ,f}\ a < c a .p\ a h^ < 

if c Qi /3 < 0. Putting now 

gl j) (X):= ]T C ^ A " and 5i j) (A):= ]T c a>/3 A a 

we have 

#(A)<^< 5 «(A) 

by construction. Moreover, again by construction, the degree of and of g^ is strictly less than 
the degree of . □ 

The above theorem tells us that the error term v x can be bounded by polynomials of strictly 
smaller maximal degree than that of 5(A). As the following example shows, the degree of gj^ and 
of g^ can in fact be much smaller than that of g^\ 

Example 3.8 Let us consider again the semi-algebraic set given by 

9 (y)=( k v l 

As can be easily checked, each integral point v in this semi-algebraic set can be written as v — g(X) 
for X 6 showing that the correction term v x is in this case. □ 



This leads us immediately to the questions of when is v x = or of when is v x e C n Z"? In both 
cases, of course, T = {g(X) : A e Z|} would be an integral function basis for C with our desired 
property T C C HZ™. 

We believe that research in this direction will make it possible to design novel algorithms for 
polynomial integer programming based on reformulation techniques. 
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